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Abstract 



Using the Cottingham formula, we give an estimate of the electromagnetic mass 
splitting of pseudoscalar heavy mesons in the beauty and charm sector. We include 
in the dispersion relation the Born term, the 1~ resonance and the positive parity 1"*" 
resonance. We also evaluate the contribution to the mass difference from the isospin 
breaking quark mass differences. Our results: — m^o = —0.83 it 0.34 MeV 

and m£)+ — m£,o = +4.33 it 0.37 MeV, are in agreement with the experimental 
measurements: 7X1^+ — m^o = —0.35 it 0.29 MeV and m£)+ — m^jo = +4.78 ± 
0.10 MeV. We also compute the mass differences in the infinite heavy quark mass 
limit, which show small deviations from the finite mass results for the B case and 
30% effects in the charm case. 



1 Introduction 



The mass difference between and mesons is an interesting physical quantity, whose 

precise Icnowledge might be of primary importance at the future B-factories. As a matter 
B (AS) > B^B^) 

of fact, the ratio ^. ^ , - , determines the relative abundance of neutral and 

BR{T{AS) B+B-) 

charged B mesons produced at such accelerators and is strongly dependent on the B^ — B^ 
mass difference, since the B pair production threshold is very close to the T(4S') mass. 

The experimental determination of 6m{B^ — B^) changed significantly during the last 
ten years, from the value 5m = —2.0 ± 1.1 ± 0.3 MeV |l[] to the values 6m = +0.9 ± 1.2 ± 
0.5 MeV (ARGUS) and 5m = +0.4±0.6±0.5 MeV (CLEO) 0. A recent measurement 
by the CLEO Collaboration gives a negative mass difference Sm = —0.41 ± 0.25 ± 
0.19 MeV, and the combined CLEO-ARGUS result is 6m = -0.35 ±0.29 MeV §. Such a 
small value has to be compared to the analogous figure for the — mass difference: 
5m{D+ - D^) = +4.78 ± 0.10 MeV. 

In the theoretical understanding of these values an important role is played by the 
isospin symmetry breaking effects related to the current quark masses. According to the 
modern picture, which incorporates the old tadpole mechanism of Coleman and Glashow 
0, the strong isospin breaking is due to the intrinsic u — d mass difference. By making the 
d quark heavier than the u quark, one can explain, at least qualitatively, all the known 
meson and baryon electromagnetic mass differences 0. In particular, the large value for 
the — mass difference can be explained by the combined effects of the u — d mass 
difference and the repulsive Coulomb energy between the c and d quark. 

According to heavy quark symmetry, the effect due to the u — d mass difference 
(qm) is independent of the heavy quark mass; therefore, in the case of B^ — B^ the 
quark mass term gives a large negative contribution and would cancel out the repulsive 
Coulomb electrostatic energy resulting in a small B^ — B^ mass difference. Such a simple 
picture, however, has to be implemented quantitatively, and this is the aim of the present 
letter. We begin by giving in Section 2 an estimate of the contribution to 6m arising 
from the u — d mass difference, using SU (3) flavour symmetry and data on Bg — B mass 
differences. Since the tiny B^ — B^ mass difference arises from the sum of two terms, 
comparable in size (a few MeV), but opposite in sign, it is desirable to have an estimate 
of the electromagnetic contribution as accurate as possible. This task is afforded by 
using the covariant Cottingham formula [§, a method employed for the calculation of 
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the electromagnetic mass differences in the hght hadron sector By the Cottingham 
approach one relates the electromagnetic (e.m.) mass difference to the forward Compton 
scattering amplitudes Ti and T2, which satisfy dispersion relations (DR) and can be put 
in a form which contains integration over space-like photon momenta = —Q^ < 0. 

The application of the Cottingham formula to the evaluation of electromagnetic mass 



differences has a long story |T2|. Previous (prior to QCD) attempts to use the Cottingham 
formula for evaluating electromagnetic hadron mass differences encountered two problems: 
the first one is the convergence of the integral and the second one is the convergence of 
the DR satisfied by Tj. The current approach to these problems involves a cut-off of the 
integral at a maximum value Q'^^x — /^^) where /i represents a scale coinciding with 
the onset of the QCD scaling behaviour this point is discussed in Section 3. As for 



the convergence of the DR, the different contributions to Im Tj can be related to different 
Feynman graphs of an effective theory including hadrons and photons. For light mesons 
this can be done by using chiral perturbation theory, and recently some determinations 
of 6m{7T^ — 7T^) and 6m{K~^ — K^) by chiral perturbation theory have appeared in the 
literature |Jl^. In this approach, the subtraction constant can be computed directly from 



the Feynman amplitudes. A similar effective theory was developed also for heavy mesons 
(for a review see |[T^ ) , and we will use it in our description of the electromagnetic coupling 



of the heavy mesons involved in the calculations (the low-lying B and B* mesons and the 
first excited positive parity resonances). Therefore, also in our approach the subtraction 
constant in the DR is directly evaluated from the Feynman amplitude. These points are 
discussed in Section 4. 

We conclude the paper by computing in Section 5 the meson mass differences in the 
rriQ —>■ 00 limit. This calculation allows a remarkable simplification of the formalism, with 
a clear view of the mechanism producing 6m. We find that the infinite limit can be well 
applied to the B case, whereas in the charm case the deviation due to the finite heavy 
quark mass is of the order of 30%. 



^Previous attempts to estimate the electromagnetic contribution to Sm{B^ — B^) used the quark 
model and QCD sum rules in |Q the result 5m{B^ — 5°) = — 1.5McV, including quark mass 
effects, was obtained. This issue has been investigated using the Cottingham formula in In [|o| 

the elastic constribution is considered with a different treatment of the light quark currents, diregarding 
inelastic contributions. In the calculations are carried out in the Nc 00 limit. The numerical 
results are in agreement with the results of this paper. 
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2 Quark mass contributions 



The contribution 6m{B^ — -B°)qm and 5m{D^ — -D°)qm from the strong isospin breaking 
u — d quark mass difference: mu — md, can be computed observing that the (approximate) 
SU{3) flavour symmetry allows us to write: 



qm 



{niu — rrid) < B^\mi\B^ > ~ (m^ — m^) < Bs\ss\Bs > 



Considering the quark mass contribution to the B^ — B^ and B^ — B^ mass differences 
we have similarly 

ITLu + rud 



5m\Bs - B+\^ + 5m\B, - 



qm 



< BAss\B, > . 



(2) 



For 6m{Bs — B) we can assume that the quark mass contribution basically coincides with 
6m, since it is of the order of the strange quark mass (~ lOOMeV), i.e. much larger than 
the expected electromagnetic mass difference (of the order aAqcn — a few MeV). Writing 



W(S, - B)qrn ~ W(B, - B) 

we obtain 

6m\B+ - 5°)q^ = \6m\Bs - B+) + W(£, - 5°) 



rriu - TRd 



2ms - (rriu + md) 



(3) 



(4) 



A similar formula holds for 5m^{D^ — D'')qm. Using experimental data for 6m'^(Bs — B^), 
6m^{Bs - 5°), 6m^{Ds - D+), 6m^{Ds - D^) § and the result given in |6l for the (scale 
independent) ratios of current quark masses from chiral perturbation theory: 



ms - {mu + md)/2 
rrid - mu 



40.8 ±3.2 



we obtain: 



5m{B+ - B° 
5m{D+ - 



/qm 



/qm 



-2.23 ±0.27 MeV 
±2.54 ±0.21 MeV 



(5) 



(6) 
(7) 



/qm 



In the limit mq — oo we expect 6m{B^ — B^ 
a prediction which is supported by the results (|^J^). 



-6m{D^ — -D°)qm regardless of mg, 
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3 The Cottingham formula 



Let us consider the mass splitting of heavy mesons due to the electromagnetic interaction. 
To be definite, we consider the B meson; its electromagnetic mass shift can be derived by 
computing: 



26 



Sm = — 



2 J (27r)4 q^ + ie ' ' 

where 

T^"'{q,p) = I J (fxe-'''^ < B{p)\T{J^{x)r{Q))\B{p) > ; (9) 

is the electromagnetic current. The Compton amplitude can be decomposed in terms 
of gauge invariant tensors: 



T'^^{q,p) = DrT,{q\ u) + Di^''Uq\ v) 



[y = V ■ q)-, where 



^2 



q^'q'' 




(10) 

(11) 
(12) 



TTiB being the meson mass. 



The first step in the calculation of the integral (P) consists of a rotation in the complex 
plane and a change of variables. Let us consider the meson rest frame, u = rriBqo- Since the 
singularities in T'^" are located just below the positive real axis and just above the negative 
real axis in the complex go plane, the integration over go may be rotated to the imaginary 
axis qo = iko without encountering any singularity. After this transformation, the integral 
involves only spacelike momenta for the photon, i.e. g^ = —Q^ = —(go + q^)- After a 
change of variables from (|q|, go) to (Q^, fco), one obtains the Cottingham formula M: 



6m 



167r3 Jo Q2 



X 



-3Ti{-Q\zko) + (l 



Q 



(13) 



In eq.(|T3D we have introduced a cut-off in the Q integration at Q^^x = • Its origin 



is as follows (see |jT3| for a detailed discussion). To take into account possible ultraviolet 
(UV) divergences, the Cottingham formula has to be renormalized. The renormalization 



is accomplished by a regularization of the integral and the inclusion of counterterms 
in the lagrangian describing electromagnetic and strong interactions of quarks, gluons 
and photons. Both the strong coupling constant and the quark masses ntg have to 
be specified at a renormalization mass scale /i; since the counterterms cancel the infinite 
contribution induced by virtual particles with momenta larger than yU, the net effect is 
analogous to a cut-off of the integral at Q'i^ax = A^^- There is a residue smooth 
dependence on /i, but it should be canceled by the /i— dependence of the renormalized 
quark masses and strong coupling constant. Typical values of /i are in the range of 1-2 
GeV, corresponding to the onset of the scaling behaviour of QCD. The presence of heavy 
quarks does not change this procedure since the relevant mass scale, in the infinite heavy 
quark mass limit, is the residual energy release and the onset of scaling is again at a few 
GeV in this variable. 

The Compton amplitudes Tj [i = 1, 2) satisfy dispersion relations (DR) in the u = p-q 
variable with Ti requiring one subtraction as follows: 



IT Jo 



dz/'2 ImTi(g2,z/'^ 

^2 



1 f'^ 
TT JO 



yl2 _ y2 



(14) 
(15) 



By employing these DR, the integral over in eq.(|I^ can be performed explicitly, with 
the result: 



+ 



a ft' 



Air Jo 



,'2 



m\Q'^ 



W2{-Q\p')k2 



-Ti(-Q^0) + 3 
2 Jo u 



/2 



/2 



(16) 



where 



1 



TT 



ImT,(g2,z/) = W,{q\u) 



(17) 



and 



My) 
My) 




(19) 
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4 B and D meson electromagnetic mass differences 



In order to evaluate the DR (|n[),(|T5[) we consider the contribution of the Born term (the 
B meson), the = 1~ resonance B*, and the positive parity resonance = 1"^ Bi. We 
notice that the Born term pole and the B* belong to the supermultiplet sf = of the 

Heavy Quark Effective Theory (HQET) (sf is the total angular momentum of the light 
degrees of freedom), whereas Bi is the = 1+ partner of the sf = supermultiplet 
(the other partner, with = 0~^, has no electromagnetic coupling to the B meson). Let 
us observe explicitly that we do not introduce the sf = (^|^ supermultiplet of heavy 
mesons containing the = and = 2+ states, for which we do not have sufficient 
phenomenological information at the moment. 

To compute the electromagnetic contribution to the B meson mass difference, we 
consider the following matrix elements {q = p' — p): 

< B{p')\J^B{p) > = f{q'){p + pr (20) 
<B*{p',e)\J!fjB{p)> = iKq')e^^^^ elQpP^ (21) 

<B,{p',e)\J^jB{p)> = -J-[K,[g^^{p''~ml)-q-{p + pr] 

+ i^2[gV'^-gV]] 4 (22) 

where e is the B* or Bi polarization vector and /, h, Ki, K2 are electromagnetic form 
factors. In general they contain two terms, describing the couplings of the electromagnetic 
current to the heavy Q = b,c and light q = u,d,s quarks, respectively: 



Co 



1 — q'^/m 



V 



h{q') = T^e(^)+ . \, (24) 

e„a 



''VJ 

K,{q') = 2eQry,{u)+ ;f (25) 

L — q I Illy 

K^iq') = 2eQri/2(^) (26) 

where uj = v ■ v', and v and v' are the heavy particle four velocities. We note explicitly 

2 1 

that, e.g. for 5+ = ub, one has Cg = -, eg = 65= +-. ^(cu) is the Isgur-Wise form 
factor 113 and ti/2{uj) is the analogous form factor describing the transitions between the 
(0~, 1^) and the (0"*", 1"^) doublets of heavy mesons p|. From HQET Jl^, at the leading 
order in l/mq: 

< B{v')\h^'h\B{v) > = mBiv^" + v'^) ^{u) (27) 
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<B*{v\e)\h^h\B{v)> = imBi{uj)e''^'"' elv'^v^ (28) 
<Bi{v\t)\h''h\B{v)> = 2^mBmB,nn{^)[{l - uj)g^'' + v'^'v'']tl . (29) 



We can write ^{u) as 



2p2 



(30) 



using the normalization condition ,^(1) = 1. The experimental determination of the 
slope contains several uncertainties (see for example the discussion in [0). A value 
= 1 ± 0.3 encompasses most of the theoretical predictions, while being in agreement 
with the data pO]. Therefore we shall take in the following = 1, which means that we 

can take for the Isgur-Wise function the following expression (with u = 1 — - — — in this 



2m| 



case 



n2 



l + UJ 



(31) 



[1 -g2/4m|]2' 

For Ti/2{uj) we take the QCD sum rule results given in we shall discuss the uncer- 
tainties related to this choice below. For the light quarks part of the electromagnetic 
current, we assume Vector Meson (p, u) Dominance of the form factor; under this hy- 
pothesis the constants A and a can be estimated as follows: Aq = niB, Ag ~ 0.5 GeV p2| . 



2V2 



9vh 



m 



V 



\fi\^/mBmB^ — 2.7 (here gy — 5.8 , fy — 0.17 GeV^, my is the p meson 



mass, and |p| ^0.1 GeV ^ parametrizes the BBiV vertex ||T5| 



Using the matrix elements and the coupling constants just introduced, we can calculate 
the electromagnetic contribution to the mass splitting of heavy mesons. The contributions 
of the different terms to the DR are as follows; the subtraction term Ti{q^, 0) is given by: 



Tl(g^O) 



-2 



where 



q^ + - m^^ 



(32) 



(33) 



As for the two structure functions ^^1,2(0'^; ^) that appear in the dispersion relations for 
Tj, they are given by: 



Wiiq',u) 



2 



hW)-hl{q') 



T 2 
rrin 

4 ^ 
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4m| 



5 - 4) 



X 



X 



l2 



+ 



4 2 

q ni Q 



^Rq 



hW) - hliq') 




g'(Xi,o + i^2,o) -2i/iiKi,o 

4 



(34) 



X 



q\K- 



1,0 



}(35) 



where /+, Kj^^ refer to S+ (resp. and and /io, /o, -^0,0 (resp. D^). 

The 1~ resonance is quite narrow (less than 1 keV); on the contrary, the 1+ axial 
vector resonance is broad enough to require the convolution of the mass difference term, 
depending upon mg^, with a lorenztian distribution centered on ms^^aver = 5.732 GeV. 
Experimental data suggest a width F = 145 MeV. The D case is computed in full analogy 
with the B one. The numerical results of this analysis are reported in Table I, for ji — 
1 GeV and (in parentheses) ji — 2 GeV. It may be useful to stress that the results are 
remarkably insensitive to variations of the cut-off // ; for example varying // in the range 
fjL — 2 — 5 GeV introduces an uncertainty of less than 8%. Another possible source of 
error is in the slope of the Isgur-Wise function p^. We find an uncertainty of ±1% for 
Sm{D'^ — D'^) and negligible for 5m{B^ — B^) when varies between 0.80 and 1.20. Also 
the uncertainties related to the choice of T1/2 are negligible, given the smallness of the 1"^ 
contribution, and we do not expect significant contributions from the = poles. 

To compare our result to the experimental data we have to add the quark mass contri- 
bution computed in Section 2; the different terms and the total theoretical prediction are 
reported in Table II which shows a good agreement with experiment within the errors. 



5 Electromagnetic mass difference in the mq oo 
limit 

We wish now to evaluate the electromagnetic mass difference — B^ in the infinite heavy 
quark mass limit, which allows a remarkable simplification of the formulae and a deeper 
understanding on the underlying physics. 
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In the nil, oo limit, since 

6m'^{B+ - 5°) = 2mB 6m{B+ - 5°) , (36) 
we get, in the m^, — > oo hmit, from previous formula, the result: 

6m{B^ - i?°)em SniBorn + Smy + SfUsuhtr {frib ^ oo) , (37) 
where ^mBorn is the contribution from the Born term and is given by 



^ amy 

OrriBorn = — 

OTT 



u u/mv 
5 arctan h 



my l + yU^/my 



(3^ 



As explained above, my — 770 MeV is the p mass and p is the cut-off; for = 1 GeV 
and n = 2 GeV, we get ^mBom = 1.5 MeV and 1.9 MeV respectively. The remaining 
contributions in (|37|) arise from the subtraction term Ti(— Q^,0) : ^mgubtr, and from the 
vector meson 1~ dispersive contribution to Wi and W2- Smv (the contribution from the 
l"*" pole vanishes). The two terms are: 

a/i^ rriB 



SttA^ 1 + p^/rriy 

2 



ajj, rriB Oim 



87rA2 1 + /iVm^ 127rA2 



H /i/my 

arctan ■ 



my l + jj?/my 



(40) 



(Aq ~ 500 MeV is the hadronic scale defined by eq. (P^)). It is interesting to observe that, 
while individually the subtraction contribution (PP]) and the vector meson contribution 
(|40|) diverge in the infinite heavy quark mass limit, their sum is finite. Therefore, summing 
up the three terms (m;, oo), we obtain: 



5m{B- 



B' 



amy 
Gtt 



5 - 



my 
2AI 



/i 



m 



V 



p/my 



arctan — h 1 + . ^ i ^ , o 

my 2A2 ; 1 + /iVmf. 



(41) 



This gives, at /i = 1 GeV and /i = 2 GeV, 6m{B+ - B^)^^ = 1.36 MeV and 1.59 MeV 
respectively, which is remarkably close to the value obtained at finite mass and reported 
in Table II. A similar formula holds for 5m(Z}+ — -D°)em in the same limit {rric — * oo): 



6m{D^ - L>°)e 



amy 



7 



my 



arctan 



my 



1 + 



p/my 



2KI) l + /iVm 



(42) 
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where the differences with ( ^ ) are only due to quark charge factors. Numerically we find, 
at fi = IGeV : 6m{D+ - D'^)^^ = 1.92 MeV (this value is 2.72 MeV for /i = 2GeV). 
These results, valid for rric —>■ oo, show significant deviations from the finite mass result 
reported in Table I. 

Besides showing the exact cancellation of the divergent term in ( |39D and (^OD, which 
confirms the scaling law 6m —>■ const {rrib —>■ oo), the previous analysis is interesting 
also because it explicitly shows the small dependence of the m?, oo results on the 
renormalization scale /x. 

We also remark that, although the 1~ state makes only a small contribution to the 
electromagnetic mass difference, its contribution seems to increase with the cut-off, as 
seen in Table I. Actually, its value at a large fi , e.g, at 2 GeV , should be smaller than 
the values we give in Table I, since the form factors /i(g^) should be further suppressed 
at large by perturbative QCD effects such that the cross sections for the production 
of 0~1~ pair in e'^e~ collisions (e.g, e~^e^ vrp) will not grow too fast with energy. This 
suppression also guarantees the convergence of the integral for the Cottingham formula 
O and make our results insensitive to the value of the cut-off /i. 



In conclusion, we can say that the small mass difference — can be understood, 
in the mq — * oo limit, as a sum of two contributions of opposite sign and similar size 
that remain finite in this limit. The electromagnetic contribution has been computed by 
the Cottingham formula and has a small dependence on the renormalization mass scale 
/i. The HQET results are very similar to those obtained at finite h mass. In the case of 
D'^ — the contributions have the same sign and add up; in this case numerical results 
show deviations of ~ 30% as compared to the predictions obtained in the HQET limit. 
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Table Captions 



Table I 

Different contributions to the electromagnetic mass differences in the B and D systems 
(units are MeV). The first value is obtained using fj, = 1 GeV ; the second value (in 
parentheses) using /x = 2 GeV. 

Table II 

Electromagnetic and quark-mass contributions to the mass differences in the B and D 
systems (units are MeV) compared to the experimental data |Q. The e.m. value is an 
average between the results obtained with /i = 1 GeV and /i = 2 GeV. 
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Table I 



5m 




Born 


1- 


1+ 


total 


6m{D+ 


-DO) 


1.72 (2.28) 


-0.09 (-0.34) 


0.004 (-0.007) 


1.63 (1.95) 


5m{B+ 




1.50 (1.87) 


-0.18 (-0.41) 


0.005 (0.01) 


1.33 (1.47) 



Table II 



5m 




e.m. 


quark mass 


total 


exp. [|] 


5m{D+ 




+1.79 ±0.16 


+2.54 + 0.21 


+4.33 + 0.37 


+4.78 + 0.10 


5m{B+ 


-fiO) 


+ 1.40 ±0.07 


-2.23 + 0.27 


-0.83 + 0.34 


-0.35 + 0.29 
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